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We present theoretical and experimental analysis of an interferometric measurement of the in-situ 
phase drop across and current flow through a rotating barrier in a toroidal Bose-Einstein condensate 
(BEC). This experiment is the atomic analog of the rf-superconducting quantum interference device 
(SQUID). The phase drop is extracted from a spiral-shaped density profile created by the spatial 
interference of the expanding toroidal BEC and a reference BEC after release from all trapping 
potentials. We characterize the interferometer when it contains a single particle, which is initially 
in a coherent superposition of a torus and reference state, as well as when it contains a many- 
body state in the mean-field approximation. The single-particle picture is sufficient to explain the 
origin of the spirals, to relate the phase-drop across the barrier to the geometry of a spiral, and to 
bound the expansion times for which the in-situ phase can be accurately determined. Mean-field 
estimates and numerical simulations show that the inter-atomic interactions shorten the expansion 
time scales compared to the single-particle case. Finally, we compare the mean-field simulations 
with our experimental data and confirm that the interferometer indeed accurately measures the 
m-situ phase drop. 


I. INTRODUCTION 

Atomtronics focuses on the creation of atomic ana¬ 
logues to electronic devices. Analogues to several elec¬ 
tronic components, such as diodes and transistors, have 
been proposed [T], while several other circuit elements 
have been experimentally realized, including capaci¬ 
tors [HE] and spin-transistors [4]. The atomic version 
of the rf-superconducting quantum interference device 
(SQUID) has been realized [SHI], and initial experiments 
towards the creation of a dc-SQUID have been performed 
[Hiin]. Both SQUID devices are formed using a toroidal 
Bose-Einstein condensate and contain one or more ro¬ 
tating weak links or barriers. Furthermore, creation of 
an atomic rf-SQUID in a ring-shaped lattice has been 
proposed mm- Theoretically persistent current states 
in (quasi) one-dimensional toroidal geometry have been 
studied extensively [MIS]. Weak links, whether super¬ 
conducting or atomic, are characterized by the relation¬ 
ship between the current through and the phase across 
the barrier |16] . Accurate measurement of this current- 
phase relationship in the atomic system is crucial for the 
characterization of atomtronic devices. 

Measurement of the in-situ phase of a condensate 
through interference is a common tool in modern cold- 
atom physics. Since the first interference between three- 
dimensional condensates was demonstrated in 1997 [m 
I18j . several experiments have used interference to infer 
details about the in-situ phase profile of condensates m- 
Vortices in condensates [ID] and fluctuations brought on 
by the two-dimensional Berezinskii-Kosterlitz-Thouless 
phase transition [2T] have also been detected interfero- 
metrically. Interference between two molecular BECs [H] 
and BECs on an atom chip have also been observed [53] . 
Recently, interference measurements have been extended 
to determine the persistent current state in a toroidal 
condensate [MJ [DS]. Reference [DS] also measured the 


current-phase relationship of a BEC in a toroidal trap 
with a rotating barrier, the atomic analogue of an rf- 
SQUID. 

In the experiment of Ref. [IS] a single condensate was 
created in a simply connected trap and subsequently split 
into two condensates. One condensate was confined in a 
toroidally shaped “science” trap and the other conden¬ 
sate was confined in a concentric disc-shaped “reference” 
trap. We will refer to these together as the “target” trap. 
A schematic and an in-situ image of atoms in a target 
trap is shown Fig. [^ The science and reference traps 
were separated by more than 5 /rm, thus atom tunneling 
between them is negligible and the condensates dephase 
rapidly because of imperfections in the splitting proce¬ 
dure. Hence, when the two condensates expand and in¬ 
terfere after turning off all trapping potentials, their rela¬ 
tive phase is random, thus representing a self-heterodyne 
measurement |26j . Rotating weak links are only applied 
to the condensate in the science trap and the other con¬ 
densate is a phase reference. The current through and 
the phase drop across the barrier were inferred from the 
spiral-shaped modulation in the density prohle for short 
expansion times. The number of spiral arms determines 
the winding number of the persistent current state, while 
their chirality determines the direction of atom flow. 

In this paper, we study in detail the interference pat¬ 
terns that result from interfering a toroidal condensate 
with a reference condensate and verify the interferomet¬ 
ric technique used in Ref. |25j to measure the current- 
phase relationship. We first study analytically and nu¬ 
merically a single-particle version of the atomic rf-SQUID 
in section [llj We find that the experimentally observed 
spirals are a short time phenomenon and both the cur¬ 
rent through and phase drop across the barrier follow 
from the geometry of the spirals. For longer expansion 
times, the spirals become modulated with concentric cir¬ 
cles due to self-interference of the torus and it becomes 


2 



FIG. 1. Panel (a) shows a schematic of atoms in a target 
trap. The inner disc and the outer ring are the reference and 
science condensates, respectively. A blue-detuned laser forms 
a rotating weak link and is shown by the blue ellipse. Panel 
(b) is an in-situ image from the experiment of atoms in a 
target trap. 


difficult to read out the in-situ phase drop. In Sec. Ill we 
describe details of our experiments with sodium conden¬ 
sates in a target trap. In addition, this section describes 
the numerical techniques used to simulate the mean-field 
Gross-Pitaevskii equation, which quantihes the effects of 
atom-atom interactions on the expanding condensates. 
Estimates of bounds on expansion times, where spirals 
can be observed in the density profile, are also derived. 
Finally, a comparison of theoretical and experimental 
results in Sec. IIVI validates the interferometric method 
for measurement of the current-phase relationship of an 
atom-SQUID. 


II. SINGLE-PARTICLE PICTURE 

We begin our study of the interference by deriving an¬ 
alytic expressions for the free expansion of a single atom 
of mass m released from a target-trap interferometer and 
give an intuitive explanation of the origin of the spirals 
in the interference pattern. To generate the interference, 
we assume that the wavefunction of our single particle 
is in a superposition of a wave localized in the reference 
and science regions, respectively. 


A. Particle in a rotating torus 

In order to solve for the wavefunctions, we first de¬ 
scribe the target trap in cylindrical coordinates x = 
{r,6,z). The science and reference traps are assumed 
to be parabolic in the radial direction, and centered at 
rs and the origin, respectively (see Fig. Ha)). The har¬ 
monic oscillator lengths are as and an, respectively. The 
common transverse confinement is harmonic with oscilla¬ 
tor length £z. We assume an « as and as^aii,£z rs- 
In addition, the science trap has a barrier or weak link 
rotating at angular frequency £1 inducing atom flow. For 
simplicity, we model the barrier in the science trap as 
a Dirac delta-function Vb{x,t) = Uow{r)d{9 — fit) with 
strength Ug, time t, and w(r) is a window function which 
is one around the radial position of the science trap and 


zero everywhere else. 

In the frame rotating with the barrier the atom is pre¬ 
pared in the time-independent state 'l'init(ai’) = ('0fl(^) + 
'ips{x))/^/2, where the iii{x) = r]i{r)(pi{9)4>ziz) are sep¬ 
arable wave functions of the science {i = S) and refer¬ 
ence {i = R) trap. Here, 4>z{z) is the unit-normalized 
ID ground state harmonic-oscillator wavefunction and 
rjiir) = jM is the radial wavefunction, 

where A/" is a normalization constant. The overlap be¬ 
tween the V’i(^) is negligible. 

The angular functions are for the refer¬ 

ence trap and the ground state of the Schrddinger equa¬ 
tion [—(P/d9‘^-\-2iKd/d9-\-U5{9)](ps{9) = E(ps{9) for the 
toroidal trap with rotating barrier. Here, k = D/Dq, 17 = 
Uq/Sq-, Do = 2£o/h, and Sq = /{2mr‘^)) « fi,^/( 2 mr|) 

is the natural energy scale of the science trap, where the 
bracket (•) indicates an expectation value over r and z 
and h is the reduced Planck’s constant. The function 
ips{9) is periodic on 9 G [— 7 r, 7 r] and a superposition of 
exp[i(/c± VT? + K^)0] with energy E = where 

e(K) is periodic in k with period one. Examples of the 
phase and magnitude of fsi9) are shown in Fig. For 
most K the phase of ^psiP) changes nearly linearly with 
9. Only for k « 1/2 and, in fact, near any half-integer 
K it changes rapidly near the barrier at 0 = 0. This 
rapid change around 9 G (—9o,9o) is accompanied by a 
decrease in density. The phase jump is 7 r(— tt) for k just 
above (below) 1/2 and the density is zero at 0 = 0 for 

K = 1 / 2 . 

We define the phase drop across the barrier as 7 = 
27r(n — s), where n is the winding number, which for the 
ground-state of the single-particle wavefunction equals to 
the integer closest to k, and the slope 


s = ^ [arg(</j(0))] 


9 — — 7T 


( 1 ) 


A graphical representation of 7 for k = 0.51 is shown in 
Fig.H In the rotating frame, the angular current 


J{k) = rsDo|</3s(6')|^ 


darg[y)5(0)] 

d9 


( 2 ) 


for any 9 and we used the fact that (1/r) = l/rg. 


B. Single-particle interference 

After turning off the target trap the atomic wavefunc¬ 
tion , '!'(£, t), freely expands and interferes. At time 
t after the release, it is imaged along the z axis lead¬ 
ing to the observable n{r, 9, t) = dz |'I'(x, t)p, where 
'^{X,t = 0) = T'init(^)- During the expansion, the 
wavefunction of the torus and disc remains separable in 
the z direction, i.e. ^i(a;, t) = Xi(r, 0, t)(/z(z, t). Thus, 
n{r,9R) = \xR{r,9,t)+xs{r,9,t)^ as / dz|(/^(z, t)p = 1 . 

It is convenient to first follow the expansion with a nu¬ 
merical solution of the Schrodinger equation in the (r, 0) 
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FIG. 2. Phase (panel a) and magnitude (panel b) of the 
single-particle ground-state angular wavefunction ips{d) as 
a function of 0 for various values of rotation rate k. The 
wavefunction is calculated in the frame rotating with a delta- 
function potential of strength Uo = 1 located at 0 = 0 . 
For K « 1/2, a sharp change in the phase occurs in the k- 
dependent region 9 € (—9o,9o). The figure also shows the 
phase-drop 7 , defined in the text, for k = 0.51. 



FIG. 3. Numerical simulation of the integrated particle 
density n(r, 9, t) of a single particle, with winding number 
equal to one, expanding in the rotating frame after release 
from a target trap. Panel (a) shows n(r, 9, t) with spirals at 
an early expansion time t = 0.25rc, evaluated at r = rs- 
Panel (b) shows a later time t = 1.25tc, where the spirals 
are superimposed with circles due to self-interference of the 
toroidal wavefunction. The density near the center has been 
truncated for better contrast. The trap parameters are an = 
0.025rs, as = O.OSrg, 17o = 1 and k = 0.51. The length of the 
sides in panels (a) and (b) corresponds to 5.12rs and 12.8rs, 
respectively. The parameters are chosen such that the overlap 
between the expanding science and reference wavefunctions is 
sufficient to show the spiral over a large range of radii. 


plane for k near 1/2. Figure|^shows n(r, 9, t) for two dif¬ 
ferent expansion times. (Time propagation was carried 
out by switching to momentum space, applying appro¬ 
priate time-dependent phase factors, and returning back 
to coordinate space.) We observe that as soon as the 
wavefunctions of the two traps overlap, the interference 
pattern consists of spirals. Later on, the self-interference 
of the science wavefunction yields circles superimposed 
on the spirals. 

We confirm these interference patterns with an asymp¬ 
totic expansion and study the associated time-scales. 
The time evolution of the reference state is 


Formally, expanding wavefunction xs(^j t) evolves as 


xs(g 



dr' r'rjs{r') 


X f de'G{r,e,r',9',t)^s{9'), 
J — TT 


(4) 


where the free-particle Green’s function m in two- 
dimensions is 


m 




( 5 ) 


exp ■ 


im [r^ + r'^ — 2rr' cos{0 — 0')] 


2ht 


where cr^(t) = + iht/m and N'i{t) normalizes the 

wavefunction. Hence, for t 3> man/H the spatial extent 
of the reference wavefunction, •\/|cr|j(t)|, is proportional 
to the expansion time, corresponding to ballistic expan¬ 
sion. In contrast, the expanding science wavefunction is 
not analytically solvable. We can, however, derive an 
asymptotic series based on the pertinent timescales of 
the expansion of the science wavefunction. The shortest 
time scale is the ballistic time tb = mag/h determined 
by the initial radial width. In addition, as will become 
clear later, there are two position-dependent timescales: 
an intermediate timescale Tc(r) = mas{r + rs)/h and a 
long timescale Ts(r) = mrrs/h. We are interested in the 
expansion time interval tb Ts{r). Figure]^ shows 

the density profile for two such times. 


We note that the integral over r' is concentrated around 
r' = rg. Consequently, the integral over O' in Eq. |^can 
be solved by noting that the phase on the right hand 
side (RHS) of Eq. [^oscillates rapidly for t ^ mrr'/h « 
Tg{r). Then, the method of steepest-descent [28] gives 
an asymptotic series for integral over O' in powers of the 
small parameter t/Tg(r). In fact, there are two stationary 
points located at 0' = 0 and O' = 0 + tt, respectively. The 
remaining integral over r' is also solved using steepest 
descent for tb based on the small parameter as/rs- 
To leading order we hnd 


Xs{r,0,t) 


'^-ir-rs)G[2al{t)] 




( 6 ) 




/iAf2it)Vr), 


where the complex, time-dependent ag{t) = cr|(l -I- 
UItb) is the square of the width of the expanding ra¬ 
dial wave-packet and l/A/2(t) is a normalization factor. 
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The wavefunction is a superposition of two expanding ID 
Gaussians centered at rs and —rg (except for the proba¬ 
bility conserving factor l/y^). The asymptotic solution 
is valid for tb t Ts{r). This excludes the region 
near the origin, where Ts(r) is small. 

It is natural to ask whether the second term in Eq. 
is important relative to the first term. Clearly, when 
\/Ws{^)\ < r rs or equivalently t < TB{r + rs)/o's = 
Tc(r) the second term is negligible. The interference of 
the first term with the reference wavefunction t) 

in Eq. leads to spirals in the density n{r, 9, t) as shown 
in Fig.]^a). Eor t > Tc{r) the second term cannot be ig¬ 
nored and interferes with the first term. It leads to circles 
in addition to the spirals as shown in Eig. [^b). An intu¬ 
itive interpretation of Tc{r) is that it corresponds to the 
time taken by signals from both antipodal points (rs, 9) 
and (rg, 9 + n) of the initial t = 0 wavefunction of the 
torus to reach the observation point (r, 9) and interfere. 
This is the self-interference of the toroidal wavefunction. 


C. Spirals 

We are now in a position to quantify the spiral 
structure for <C t <C tc- We write Xi(r, 0,t) = 

Hi (r, 9, t) exp[z^i(r, d, t)], where ni{r,9,t) is the prob¬ 
ability density and ^i{r, 9, t) is the phase. The integrated 
density becomes 

n{r, 9) = ns{r, 9)+nR{r, 9)+2y/ns{r, 9)nR{r, 9) cos^(r, 9), 

where ^(r, 9) = &) — and we suppress the time 

argument for notational simplicity. The last term on the 
RHS of this equation describes the interference of the 
wavefunctions in the two traps. 

For the above time interval the second term in Eq. 
can be ignored, so that nR{r,9,t) is independent of 9, 
ns{r, 9, t) is a separable function of r and 9, and ^(r, 9) « 
arg[(^s(0)] — hrrs/{mt). (The argument arg[(/)s(d)] is 
defined as a monotonic function of 9.) Then, spirals 
correspond to curves of constant phase ^ (r, 9) in the 
(r, d) plane. The densities ni{r,9) only lead to a slowly- 
varying envelope in r and suppression of the signal near 
9 = 0 that is most pronounced for half-integer k. Con¬ 
sequently, a spiral is described by the parametric curve 
r(u) = (^o+ai'g['^s(^(u))])xIit/(mrs) and 0 (u) = -tt+u 
mod 27r, where is a constant (typically chosen such 
that n{r, 9) is a local extremum) and u is the free pa¬ 
rameter. In the absence of a rotating barrier but for 
a non-zero winding number n of the toroidal state, we 
find arg[(/?s(d)] = n9 and the interference pattern has 
Archimedean spirals with r{u) = (^o + nu) x ht/{mrs) 
and 9{u) = —tt -|- u mod 27r. These smooth spirals have 
been observed experimentally [Ml |M] . 

A schematic of a spiral is shown in Fig. [^ at a sin¬ 
gle expansion time t for k slightly greater than 1/2, a 
case where (ps{0) has a sharp phase jump across the bar¬ 
rier near 9 = 0. For |0| > 9q the spirals smoothly wind 



FIG. 4. Schematic of a spiral-like contour (solid line) in the 
integrated density for k slightly larger than 1/2, so that the 
winding nnmber n = 1. The contour has a constant phase 
i{r,9) = ^ 0 - The phase of <ps(^) varies rapidly in the wedge 
0 £ (—6 ^o,^o)- In addition, a Archimedean spiral (dashed 
line) with the same initial angnlar velocity as the solid line 
is shown. Its parameters as well as the lengths 5 and A are 
dehned in the text. 

around the origin. In contrast, for 9 G (—9o,9(j) there is 
a sharp, nearly discontinuous change in the spirals. For 
K away from half-integer values the spirals are smooth 
everywhere. The geometry of a spiral is completely de¬ 
termined by the phase ^(r, 9) where the number of spiral 
arms is the winding number n. The densities nR(r,9) 
and ns{r^9) determine how many windings of a spiral 
are visible along the radial direction. 

We characterize the discontinuity or jump of the spirals 
by lengths <5 and A shown in Fig. [^ The quantity 5 = 
2TTht/(mrs) is the radial fringe spacing and measures the 
increment in r as ^(r, 9) is increased by 27r at a fixed 9. 
Moreover, A = rA{u -I- 27r) — rA{u) = s x 2Trht/{mrs), 
where we used the Archimedean spiral rA{u) = {^o+su) x 
ht/(mrs) and 9a{u) = —tt -|- u mod 27r, and s is defined 
by Eq. 0 Intuitively, A is the radial distance covered by 
a spiral when it is smoothly continued across the barrier 
region. The two lengths depend on the dimensions of the 
torus and expansion time t. 

The ratio A/6 = s is independent of the radial wave- 
function and expansion time. In fact, we can interpret 
A/i5 as a measurement of the phase across the barrier 7 , 
since 

7 = 27r(n — A/(5). (7) 

Moreover, it is a measurement of angular current J{k), 
as the hydrodynamic equation Eq. j^at 9 = —tt gives 

J(k) = rsDo|(/3s(-7r)p (A/5 -k) . ( 8 ) 

For t > Tc radial rings will get superimposed on the 
spirals due to the self-interference, making extraction of 
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curves of constant ^(r, 0) more difRcult. Moreover, when 
t ^ Ts(r), the derivatives of the initial angular wave- 
function become important; finally, for t 3> Ts(r), the 
probability distribution resembles the Fourier transform 
of the initial wavefunction, which has no spirals and the 
in-situ phase can not be read out. 


III. EXPERIMENTAL ATOM SQUID AND 
MEAN-FIELD SIMULATION 


We have also performed interference experiments with 
quantum-degenerate Sodium atoms in a target trap 
as well as simulations based on the mean-field Gross- 
Pitaevskii equation (GPE). These results can be com¬ 
pared to our single-particle analysis and show the role of 
atom-atom interactions present in ultra-cold atomtronic 
experiments. 

The experimental setup is described in Sec. Ill A De¬ 


tails of our numerical methods to simulate the GPE are 
given in Sec. |IIIB[ while Sec. |IIIG| describes expansion 
timescales based on a self-similar expansion of a BEG 
from a target trap [21]. Section IV compares our results 


and enables us to verify the extraction technique used in 
Ref. |25| for the phase-drop across the barrier in terms of 
a measurement of A/i5. 


A. Experimental setup 

We have performed interference experiments in a tar¬ 
get trap. We create a ^^Na BEG in a target trap with 
approximately 7 x 10® atoms and a chemical potential, 
fj./h Ki 2^ X (2 kHz). Details of the creation of the trap¬ 
ping potential can be found in Refs. The target 

trap has an external toroid with a radius of 22.4(4) ^m 
and radial trapping frequency of 240 Hz. Its central disc 
has a flat-bottomed potential and contains about 25% of 
the total atoms. The transverse trapping frequency of 
both traps is « 600 Hz. This leads to a Bose condensate 
with a measured Thomas-Fermi radial width of about 6 
/im in the toroid and a Thomas-Fermi radius of about 

5 /rm in the disc. The barrier potential has a Gaussian 
profile with a height less than the chemical potential of 
the atoms in the science trap. Its 1/e^ full width is « 

6 /im. Persistent current states are created by adiabati- 
cally ramping up of the height of the barrier with a fixed 
rotation rate. The atom cloud is imaged along the trans¬ 
verse direction by absorption imaging, which measures 
the intensity of resonant light transmitted through the 
expanding gas. 


of a BEG with a stationary weak link or barrier but oth¬ 
erwise the same trapping potentials and atom number 
as in the experiment. We use imaginary-time propaga¬ 
tion and a two-dimensional effective Lagrangian Varia¬ 
tional Method (2D LVM) |3D1 IM] . assuming a scatter¬ 
ing length a = 2.8 nm. The method is a variational 
technique whose trial wave function is the product of 
an arbitrary function in the (r, 9) plane and a Gaussian 
in the z direction with an (imaginary-)time-dependent 
width and a phase that is quadratic in z. This Ansatz 
leads to (a) a two-dimensional effective GPE whose non¬ 
linear coefficient contains the width of the Gaussian and 
(b) an evolution equation for the width that depends 
on the spatial integral of the fourth power of the ab¬ 
solute value of the solution of the effective GPE. We de¬ 
note this solution by 'l'stat(^) and normalize such that 
/ d®i'|4'stat(^)P = N, the total atom number. In partic¬ 
ular, we can find the angular density profile of the science 
trap pstat(6*) = J' rdrdz\'i>sts.t{r,6,z)\'^, where the radial 
integral only encompasses the science or toroidal trap. 

The second step is to add the rotation of the bar¬ 
rier by multiplying the stationary (and positive) 'kstat(^) 
with a spatially dependent phase that leaves the density 
profile unchanged, i.e. 'I'Gp(a^) = 'I'stat(5’)e*^^^^. The 
phase profile C(^) is zero around and inside the cen¬ 
tral disk and near the torus only depends on 9. For a 
given rotation rate k and winding number n it is found 
by simultaneously solving the hydrodynamic expression 
J = rs^oPStRt{S)idC{9)/d9-K) and C(7r)-C(-7r) = 2mr. 
(Gompare to Eq. as well as see the supplemental ma¬ 
terial in Ref. (ISj). The solution is similar in behavior to 
those shown in Fig. and the phase drop follows from 
7 = 27r(n — s), where s = d({9)/d9 \e=-T^. 

This phase-imprinting procedure is valid as long as the 
height of the barrier is l ess than the chemical potential, 
the healing length £ = ^Jhj (2m/i) « 0.5 /im is small com¬ 
pared to the width of the barrier (« 6 /tm), and the 
speed of the barrier is small compared to the speed of 
sound c = \Jp/m. These conditions are also met in the 
experiment. 

Finally, we simulate the expansion of our BEG wave- 
function released from a target-trap by solving the (real) 
time-dependent Gross-Pitaevskii equation using the same 
2D-LVM method. The GPE solutions have only been 
modified to include the effects of absorption imaging. 
The non-zero point-spread-function of the imaging sys¬ 
tem is taken into account by convolving the simulated 
transmission with an Airy disk of the appropriate size. 


C. Expansion time scales 


B. Numerical simulation 

The initial wavefunction, 4'Gp(2f), of the condensate in 
the target trap is found in a two-step process. We first 
solve the Gross-Pitaevskii equation for the wavefunction 


References [221133 showed that a harmonically trapped 
and interacting Bose condensate expands at a much faster 
rate than an non-interacting gas of the same size. Here, 
we perform a similar analysis for expansion from a target 
trap. In fact, under the assumptions valid for phase im¬ 
printing in Sec. HIB it is sufficient to study expansion 
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FIG. 5. (Color online) Panel a) False color images of the interference pattern in the density profile after a 17 ms expansion 
time for four rotation rates of the barrier. The atom number density increases from blue to red with blue corresponding to 
zero density. Top and bottom rows show images from our experiment and GPE simulations with the same trapping potentials 
and atom number, respectively. The extracted A/d for each rotation rate is shown above the images. The winding number is 
zero for all images. Panel b) Images of experimental (top) and GPE (bottom) density prohles for four expansion times of a 
non-rotating condensate released from a target trap without a barrier. 


from a BEG in a toroidal trap without a barrier or ro¬ 
tation. We assume that the interactions are sufficiently 
strong that the Thomas-Fermi approximation holds along 
the r and z directions. The BEG wavefunction is then 
independent of 9 and the harmonic confinement in the 
toroidal trap along the r and z directions leads to a BEG 
with Thomas-Fermi radius, cttf, such that cttf ^ fs- 
Here, for simplicity we assume the same trap frequency 
along the two directions, i.e. u>r = = w. 

Immediately, after the release of the toroidal trap the 
BEG expands rapidly in the r and z directions as the 
interaction energy gets converted to kinetic energy. This 
defines a ballistic timescale tb (We use tilde to de¬ 
note timescales associated with expansion of the inter¬ 
acting BEG.) As (Ttf ^ rs, we can locally approxi¬ 
mate an angular section of the torus as a two dimen¬ 
sional tube, which expands along its transverse direc¬ 
tions. Such an elongated BEG undergoes a self-similar 
expansion [29l|32]. That is, in the hydrodynamic picture 
of the BEG and cylindrical coordinates, the density is 
n(r, z, t) « n(rs + (r — rs)/\{t), z/\(t), t = 0) while the 
velocity field vix.t) = (vr{r,t),Q,Vz{z,t)) with Vrir.t) = 
(1 - A(t)“^)(r - rs)/t and Vz {z,t) = (1 - X{t)~‘^)z/t. 
The scaling factor \{t) = y/l+UH?, which implies 
Tb = = mag/h and is the same as the single-particle 

ballistic time tb , even though the radial size of the BEG 
wavefunction utf S> as- 

For Tb, the interaction energy has been converted 
to kinetic energy, the density profile has spirals, but the 
cloud is expanding more rapidly than the single-particle 
case. Hence, we expect that the time scale, fc(r), where 
the spirals become modulated with circles due to the self¬ 
interference of the toroidal BEG, will be shorter than 
the equivalent single-particle time scale, Tc(r). We can 
derive tc following the intuitive understanding of signals 
from antipodal points (rs,9) and (rs,9 -|- tt) at t = 0 
reaching (r, 9) at t = tq. In other words, we require that 
the radial size of the toroidal BEG, A(fc)o'TF, is larger 


or equal to the distance between the observation point 
and the antipodal points, i.e. r + rs and r — rg- Hence, 
fc « (j’+'cs)/(a;(TTF) = (o's/crTF)'rC) which is smaller 
than Tc- 


IV. COMPARISON OF THE EXPERIMENT 
WITH THEORY 

We compare our experimental data and GPE simula¬ 
tions in Fig. by showing the dependence of the inter¬ 
ference pattern on the rotation rate of the barrier and 
the expansion time. Figure [^) shows typical expanded 
clouds at 17 ms expansion time from our experiment and 
simulated GPE expansions for various rotation rates of 
the barrier leading to condensates with winding number 
n = 0. Firstly, we see radial interference fringes at fixed 
9 and azimuthal interference fringes at fixed r similar to 
those in Fig. The ratio A/5 from these experimental 
images is extracted following the procedure explained in 
Fig. a The phase-drop across and the current through 
the barrier then follows from Eqs. and respectively. 
Near 9 — 0, where the barrier is located before release, 
the density profile has radial stripes, which are absent 
from the single-particle simulations and a consequence of 
interaction-induced expansion of atoms into the density 
depleted weak-link region. Lastly, star-like structures, 
which are due to residual azimuthal asymmetries in the 
toroidal potential, are visible. 

Figure [^) shows expanding, rotationless clouds re¬ 
leased from a trap without barrier for various expansion 
times. For observation radii r > 60 fim and small expan¬ 
sion times t <20 ms, the experimental data and GPE re¬ 
sults show no evidence of self-interference of the toroidal 
BEG consistent with t < fc{r)- For longer expansion 
times we observe self-interference. It is prominent near 
the cloud center, where radial fringes emerge with half 
the spacing of those at large radius. 
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FIG. 6 . Radial fringe spacing, S, of the interference pat¬ 
tern as a function of time elapsed after the release of the 
target trap. The data is for a toroidal trap without a barrier 
and a BEG without winding. The experimental, GPE and 
single-particle fringe spacings are shown by red dots with one- 
standard deviation statistical error bars, blue markers and a 
black line, respectively. The value of rs has a uncertainty, 
which is shown by the shaded region around the black line. 

In Fig.[^the size, shape and interference pattern of the 
clouds in the GPE simulations agree well with those of 
the experiment. The agreement is made quantitative in 
Fig.§ for the target trap without a barrier and a BEG 
without winding (n = 0). The hgure shows the radial 
fringe spacing, <5, from the experimental data, GPE simu¬ 
lations and the single-particle expression 6 = 2Trht/{mrs) 
as functions of expansion time. The three cases are in ex¬ 
cellent agreement, indicating that this fringe spacing is 
determined by the geometry of the system, i.e. the radius 
of the torus. 

Figure]^ shows the extracted A/S as a function of the 
imprinted phase drop 7 across the barrier for the GPE 
simulations in Fig. [^). The result agrees within our 
uncertainties with the single-particle prediction, which 
indicates that interactions do not change the phase drop 
over the barrier region even though the angular density 
profile is distorted during the expansion. In other words, 
an extraction of the phase drop from a measurement of 
A/6 is valid even when the GPE and experiment have 
radial stripes for small 6 near the weak link. The latter 
are absent from the single-particle interference pattern. 

V. CONCLUSION 

We have experimentally and theoretically investigated 
an interferometric measurement of the phase drop in an 


atomic-SQUID. The atomic-SQUID consists of a BEG in 
a toroidal trap with a rotating barrier. The phase drop 
across the barrier is measured by interference with a ref¬ 
erence disc BEG after release from the trapping poten- 



FIG. 7. The ratio A/ii as a function of in-situ phase-drop 
7 across the rotating barrier from the GPE simulations of 
Eig. [^) (markers) and the single-particle prediction (solid 
line). Error bars are one standard deviations uncertainties 
from the fit to the density profile. 


tials. We have studied the single-particle case and find 
that the structure of the interference pattern depends on 
the expansion time after release. For short times, it con¬ 
sists of spirals, which have the same number of arms as 
the winding number of the toroidal wavefunction. The 
phase along a spiral is the same as the in-situ phase of 
the angular wavefunction. Moreover, we find that the 
phase drop across the barrier and the current through it 
determine the geometry of spirals. For longer times the 
spirals get superimposed by circles making phase readout 
difficult. 

The conclusions from the single-particle model are 
confirmed by experiments with Bose condensed sodium 
atoms and numerical simulations based on the Gross- 
Pitaevskii equation even though inter-atomic interactions 
speed up the expansion, thereby shortening the associ¬ 
ated time scales. In particular, one feature that is not 
changed is the fringe spacings of the interference pattern. 

Most importantly, we have confirmed that the phase- 
drop across the barrier as measured by our experiment 
agree with those of our single-particle model and mean- 
field simulations and accurately reflect the in-situ value. 
This confirmation opens up the possibility of using this 
technique for measuring the current-phase relationship 
of, for example, excitations or weak links in degenerate, 
superfluid Fermi gases. 


[1] B. T. Seaman, M. Kramer, D. Z. Anderson, and M. J. [2] J. G. Lee, B. J. Mcllvain, G. J. Lobb, and I. Hill, W. T., 
Holland, Phys. Rev. A 75, 023615 (2007) 







Sci. Rep. 3, 1034 (2013) 

[3] S. Krinner, D. Stadler, D. Husmann, J.-P. Brantut, and 
T. Esslinger, Nature 517, 64 (2015) 

[4] M. C. Beeler, R. A. Williams, K. Jimenez-Garcia, L. J. 
LeBlanc, A. R. Perry, and I. B. Spielman, Nature 498, 
201 (2013) 

[5] A. Ramanathan, K. C. Wright, S. R. Muniz, M. Ze- 
lan, W. T. Hill, C. J. Lobb, K. Helmerson, W. D. 
Phillips, and G. K. Campbell, Phys. Rev. Lett. 106, 
130401 (2011) 

[6] K. C. Wright, R. B. Blakestad, C. J. Lobb, W. D. 
Phillips, and G. K. Campbell, Phys. Rev. Lett. 110, 
025302 (2013) 

[7] S. Eckel, J. G. Lee, F. Jendrzejevrski, N. Murray, C. W. 
Clark, C. J. Lobb, W. D. Phillips, M. Edwards, and 
G. K. Campbell, Nature 506, 200 (2014) 

[8] C. Ryu, P. W. Blackburn, A. A. Blinova, and M. G. 
Boshier, Phys. Rev. Lett. Ill, 205301 (2013) 

[9] F. Jendrzejewski, S. Eckel, N. Murray, C. Lanier, M. Ed¬ 
wards, C. J. Lobb, and G. K. Campbell, Phys. Rev. Lett. 
113, 045305 (2014) 

[10] L. Amico, D. Aghamalyan, F. Auksztol, H. Crepaz, 
R. Dumke, and L. C. Kwek, Sci. Rep. 4, 4298 (2014) 

[11] D. Aghamalyan, M. Cominotti, M. Rizzi, D. Rossini, 
F. Hekking, A. Minguzzi, L.-C. Kwek, and L. Amico, 
New J. Phys. 17, 045023 (2015) 

[12] V. A. Kashurnikov, A. 1. Podlivaev, N. V. Prokof’ev, and 
B. V. Svistunov, Phys. Rev. B 53, 13091 (1996) 

[13] Y. Kagan, N. V. Prokofev, and B. V. Svistunov, Phys. 
Rev. A 61, 045601 (2000) 

[14] H. P. Biichler, V. B. Geshkenbein, and G. Blatter, Phys. 
Rev. Lett. 87, 100403 (2001) 

[15] M. Cominotti, D. Rossini, M. Rizzi, F. Hekking, and 
A. Minguzzi, Phys. Rev. Lett. 113, 025301 (2014) 

[16] K. Likharev, Rev. Mod. Phys. 51, 101 (1979) 

[17] M. R. Andrews, C. G. Townsend, H.-J. Miesner, D. S. 
Durfree, D. M. Kurn, and W. Ketterle, Science 275, 637 


(1997) 

[18] Y. Castin and J. Dalibard, Phys. Rev. A 55, 4330 (1997) 

[19] J. E. Simsarian, J. Denschlag, M. Edwards, C. W. Clark, 

L. Deng, E. W. Hagley, K. Helmerson, S. L. Rolston, and 
W. D. Phillips, Phys. Rev. Lett. 85, 2040 (2000) 

[20] S. Inouye, S. Gupta, T. Rosenband, A. P. Chikkatur, 
A. Gorlitz, T. L. Gustavson, A. E. Leanhardt, D. E. 
Pritchard, and W. Ketterle, Phys. Rev. Lett. 87, 080402 
( 2001 ) 

[21] Z. Hadzibabic, P. Kruger, M. Cheneau, B. Battelier, and 
J. Dalibard, Nature 441, 1118 (2006) 

[22] C. Kohstall, S. Riedl, E. R. Sanchez Guajardo, L. A. 
Sidorenkov, J. Hecker Denschlag, and R. Grimm, New 
J. Phys. 13, 065027 (2011) 

]23] Y. Shin, C. Sanner, G.-B. Jo, T. A. Pasquini, M. Saba, 
W. Ketterle, D. E. Pritchard, M. Vengalattore, and 

M. Prentiss, Phys. Rev. A 72, 021604 (2005) 

[24] L. Gorman, L. Chomaz, T. Bienaime, R. Desbuquois, 
C. Weitenberg, S. Nascimbene, J. Dalibard, and 
J. Beugnon, Phys. Rev. Lett. 113, 135302 (2014) 

[25] S. Eckel, F. Jendrzejewski, A. Kumar, C. Lobb, and 
G. Campbell, Phys. Rev. X 4, 031052 (2014) 

[26] L. Richter, H. Mandelberg, M. Kruger, and P. McGrath, 
IEEE J. Quantum Electron. 22, 2070 (1986) 

[27] R. P. Feynman, A. R. Hibbs, and D. F. Styer, Quantum 
Mechanics and Path Integrals: Emended Edition (Dover 
Publications, 2010). 

[28] C. M. Bender and S. A. Orszag, Advanced Mathematical 
Methods for Scientists and Engineers (McGraw-Hill, New 
York, 1978). 

[29] Y. Castin and R. Dum, Phys. Rev. Lett. 77, 5315 (1996) 

[30] M. Edwards, M. Krygier, H. Seddiqi, B. Benton, and 
C. W. Clark, Phys. Rev. E 86 , 056710 (2012) 

[31] V. M. Perez-Garcia, H. Michinel, J. 1. Cirac, M. Lewen- 
stein, and P. Zoller, Phys. Rev. A 56, 1424 (1997) 

[32] Y. Kagan, E. L. Surkov, and G. V. Shlyapnikov, Phys. 
Rev. A 54, R1753 (1996) 







